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Abstract 



It has been proposed to determine the CKM matrix element \V u i\ in a model-independent way 
from a combination of rare and semileptonic B and D decays near the zero recoil point. An essential 
ingredient in such a determination is a heavy quark symmetry relation connecting the form-factors 
appearing in B — » K*e + e~ to semileptonic form factors relevant for B — > pev. We estimate the 
leading corrections to this symmetry relation, of order a s {m^) and A/m^, pointing out that they 
can be as large as 20%, depending on the value of the matrix element of a dimension-4 operator. 
Dimensional analysis estimates of this matrix element give a corresponding uncertainty in \V U \\ of 
the order of a few percent. 
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The Cabibbo-Kobayashi-Maskawa matrix element V ub is an important ingredient of the 
flavor structure of the Standard Model. Its magnitude determines one of the sides of the 
unitarity triangle. Although at present it is one of most poorly known CKM parameters, 
several promising methods have been proposed to extract it from experimental data (see |IJ 
for a recent review). 

We focus in this paper on a particular model- independent method for determining \V u b\ 
from exclusive rare and semileptonic B decays proposed in 0, [3|, |J]. There are two basic 
ingredients going into such a determination: 

1. Heavy quark symmetry relations between heavy- light form factors |6| connect the 
rate for B — > K*e + e~ at the zero recoil point g^ ax to the rate for B — > peu, up to 
A/m b and SU(3) breaking corrections. 



dT(B -> peu)/dq 2 \V ub \ 2 Stt 2 I \f B ^(y)\ 2 I q) 



dT(B - K*£+£-)/dq2 \V tb V t * s \ 2 a 2 \C 9 \ 2 + |C 10 | 2 \f B ^ K '{y)\ 2 1 + A(y) 

In this relation the ratio of decay rates is taken at a common value of the parameter 
y defined by v ■ pv = mvy- The correction A(y) parametrizes the contribution of the 
magnetic penguin operator Q 7 near y — 1. 

2. The SU(3) breaking corrections introduced in step 1 can be eliminated with the help 
of semileptonic decay D — ► (p, K*)ev data, by using an approximate equality to unity 
of the double ratio |7j 

First evidence for the mode B —> K*£ + £~ has been recently reported by the BaBar 
collaboration, and an upper limit was given by BELLE [§. ^| 

B(B -> K*£+£-) = (1.68+°;1 ±0.28) x (BaBar) (3) 

< 14 x 10' 7 (90% CL) (BELLE) . 

This suggests that a determination of \V u b\ using these decays might become feasible in a 
not too distant future. 

There are several sources of theoretical uncertainties connected with such an approach. 
The dominant theoretical uncertainty is connected with long-distance contributions com- 
ing from four-quark operators in the weak Hamiltonian Q x = (s Q ,7 /i Pc,6 i a)(c i g7 /i Pi / c a ), 
Qi — (sjfiPLfyicjpPLc) [TCfl . Their effect can be absorbed into a redefinition of the Wilson 



coefficient C$ ]T0| but this modification is generally process-dependent. It can be com- 
puted reliably in perturbation theory as long as the e~e + invariant mass is sufficiently low. 
Such computations have been performed for inclusive B — > X s e + e~ [[II], [TTJ and exclusive 
B — > K*e + e~ decays |T4] , |T5|] . At small recoil, as considered in this Letter, this method 
is not applicable and one has to resort instead to a phenomenological parameterization of 



these effects in terms of sums over J/ip resonances | 16| . While this is not a controlled ap- 
proximation, the validity of such an expansion can be tested by measuring other observables 
such as q 2 spectra and/or angular asymmetries j|, |i~5| . 
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A better understood source of uncertainty is introduced in step 2 as corrections to the 

double ratio R(y). From power counting, such corrections are of order ( — — ) ~ 7%. 

A more precise estimate has been performed in H using the chiral perturbation theory for 



vector mesons [ 1y| . The corresponding deviation of the ratio R(l) from unity was found to 
be small, under 1%. 

Finally, another source of corrections is introduced in the first step of the method. As 
shown in ||, the contribution of the electromagnetic penguin operator Qj to the B — > 
K*e + e~ decay rate (parameterized by A(y)) can be computed in a model-independent way 
at the zero recoil point y = 1 using heavy quark symmetry The corrections to this prediction 
are of order 0(AQCD/fnb) and were thus thought to be negligible. 

Recently, the structure of the leading corrections to the heavy quark symmetry relations 
for heavy-light decays ||, |6| has been studied in They come from dimension 4 operators 
in the matching of the weak current onto heavy quark effective theory (HQET) operators, 
and from hard gluon effects appearing as Wilson coefficients in the HQET. We are therefore 
now in a position to study the leading corrections in step 1 of the \V u b\ determination. 
We find that there are calculable corrections to the heavy quark symmetry relations of the 
order of 10%, plus an additional correction which depends on an unknown form factor of a 
dimension-4 operator. Estimates of this form factor from dimensional analysis result in an 
uncertainty on \V u b\ of the order of 3%. 

The amplitude for B — > K*e + e~ receives in general contributions from all operators in 
the weak radiative b —+ s£ + £~ effective Hamiltonian |TIJ. The result can be expressed in 
terms of the matrix elements of the three operators 

Qr = ^(m^P R b)F^ (4) 

Q9 = ^(s^P L b)(e^e) (5) 

Qio = ^(s^PLb)(e^ l5 e) , (6) 

where the effects of the remaining operators Qi-q are understood to be absorbed into redef- 
initions of the Wilson coefficients Cj s 9 10 . 

The amplitude for the decays considered here requires the knowledge of the form-factors 
of the vector, axial and tensor currents sTb. We define them as |S], [TS| (with the convention 
e 0123 = 1) 

(V(p',e)\q lfl b\B(p)) = ig{q 2 )e, vX(7 tl{p + p') x {p~p\ (7) 
(V(p',e)\q W b\B(p)) = f(q 2 K+ a + (q 2 )(e* -p)(p + p'), (8) 



-{q 2 W ■ p){p - p') 



and 



{V(p\e)\qa^ v b\B(p)) = g+(q 2 )e^ x<7 el(p + p') a + g-(q 2 )e^ x<7 e* x (p - p') a (9) 

+ h(q 2 )£^ Xa (p + p') x (p - p') a (e* ■ p) . 

The differential decay rate for B —* K*e + e~ can be written as a sum of contributions 
corresponding to well-defined helicities A of the vector meson 

±T{B -> K*e+e~) oc £ {{H^ + |^ A) | 2 } . (10) 



A=±1,0 
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The two helicity amplitudes H^' A ' correspond to the vector and axial coupling to the leptons, 
respectively. Expressed in terms of the form factors defined in (0)-©, they are given 
explicitly by 

fl£V) = C 7 2 -^( g+ ( q ")[±2rn B \q\ +m l- m l] +g 4 q 2 ) q 2 ) (11) 
+ C 9 [±2m B g(q 2 )\q\ + f(q 2 )] 

<V) = C ^^T V {9 + {<f){< + 3< - q 2 ) + g-(Q 2 )« -m\- q 2 ) + 4% 2 )m|g 2 ) 
+ C,—^ (f(q 2 )km 2 B - m 2 v - q 2 ) + 2m 2 B q 2 a + (q 2 )) (12) 



and 



tfilV) = C w (±2m B g(q 2 )\q\ + f(q 2 )) (13) 
Hi A \q 2 ) = C 10 —^-= (f{q 2 )\{m B - m 2 v - q 2 ) + 2m 2 B q 2 a + (q 2 )) (14) 

At the zero recoil point y — 1 (corresponding to maximal q 2 , g^ ax = (m B — my) 2 ), the form 
of the helicity amplitudes simplifies drastically and can be written only in terms of one axial 
form factor f(y = 1) as 

Hg\y = 1) = Hj, v \y = 1) = C fl /(1)(1 + 5(1)) (15) 
HH\y = l) = H ( A \y = l) = C 10 f(l) (16) 



where we defined 



with 



%) = ;^WcT 7m (17) 



F(y) = 9+(y)( m B + m v ) + g_(y)(m B - m v ) . (18) 

Inserting the expressions for the helicity amplitudes at the end-point (|l|), (0) into the 
rate for B — > K*£ + £~, gives the following result for the correction A(y) appearing in ([!]) 

- w^f + ©} ■ (19 > 

The form-factor ratio JF(1)//(1) is predicted from heavy quark symmetry to be 1 at 
leading order in the heavy quark limit 0, [f| . The purpose of this Letter is to estimate the 
leading corrections to this result, and study their effect on the \V u b\ determination. 

Using the results of [19| one finds the following prediction from heavy quark symmetry 

HV) = f«B + — —) f(v) + ^m B m v (l - -) g{y) + — V x {y) (20) 
\ m B J \ K\J m B 

+ 2m B m v ( 1 - A + mV ) (y - l)g{y) + 2m v (y - l)V(y) + ■■■ , 
V m B 



where the ellipses denote contributions suppressed by A 2 /ml relative to the leading term. 

The subleading form factors T>(y) and T>i(y) appearing in fl20|) are defined by matrix 
elements of the dimension-4 currents 

(V(p',e)\qiD^h v \B{v)} = V(y)ie ilvX(J El Vx v , a (21) 

iy[ji,e)\qiD^h v \B[v)) = p l£ ; + p + ( £ *.p)( P/ ,+p;) + p_(£*.p)(p M -p;). (22) 

The equation of motion for the heavy quark field iv ■ Dh v = implies a relation among the 
form factors of the qiD^^h v current, such that only two of them are independent. 

The coefficients K\ and K5 contain hard gluon corrections. The first coefficient K\ = 
—Co(mb)/c' Q (mb) is defined as the ratio of two Wilson coefficients appearing in the matching 
for = q-f^b and J'^ = qia^b 

4 ] = c ° (/^A + $ W9"A + 0{l/m b ) . (23) 

The values of the Wilson coefficients can be extracted from the next-to-leading computation 
of [2(| and their explicit values can be found in ||19|| . For our estimate we only require their 



expressions to one- loop order, which give Ki = 1 + 0(a 2 (m b )). 

The coefficient k 5 is defined analogously in terms of the Wilson coefficients appearing in 
the matching of the currents J 5m = qia ^ v v h '756 and J' 5fl = (g^ v — v^v^qj^j^b 

A'l = % (j^Tfiilshv + ^ (j^gv^sK + 0(1/ m b ) . (24) 
The explicit results for the Wilson coefficients depend on the 75 definition and are given by 

coK) = ci(m») = 1 - 4a l (mfe) + 0(a 2 s ) (NDR,HV) (25) 

07T 



, 1 + 0(a 2 s ) (HV) 
c K) = c x {m b ) = 1 [1 _ 4^ + 0(a , } (NDR) (20) 

This gives for the coefficient k 5 



_ co{m b ) _ \ l + 0{a 2 s ) (NDR) 

So(m b ) [l-^ + Oial) (HV) ^ 

For completeness we give also the corrected HQET symmetry relations used in deriving 

m m 

- 9-) + ^rn B g = -2{m v y - K)g{y) - —f{y) - 2V{y) + 0(m^ /2 ) (28) 

m B 

g + + g_ - 2m v yg - k 5 /= (29) 

m B 

-2^(Ay - m v )g{y) + —f + \{m v m B yV{y) + V x {y)) + 0{mf /2 ) 
m B m B m B 

From ( p0|) one finds for the ratio of form factors appearing in 5(1) at the zero recoil point 
y = l 

JF(1) A -m y / 1 \ o(l) 2 VAl) . . 

= K5 + YL + 2mBmv (l-—)^. + £1 30 

f(l) m B \ KiJ /(l) m B f(l) 

1-0.1 + ^--^ (NDR) 



™B /(I) 

1-0.09 -0.1 + -2--^$ (HV) 

m B /(l) V / 



This is the main result of this Letter. We turn now to a discussion of the numerical impact 
of these corrections. 

In the NDR scheme, the radiative corrections from start at two-loop order and thus 
are negligibly small (the numerical value for k^ V ^ uses a s (m b ) = 0.22). The subleading 
l/m& correction from the second term is a negative ~ 10% contribution, and is kinematically 
enhanced by the large ratio (mp — A)/mg ~ 0.10, where we took A = 350 MeV. The last 
term depends on the subleading form factor T>i(l), which has not been calculated yet. From 
dimensional analysis its contribution is expected to be of order A/m b — 0.10. Therefore, 
if this term turns out to be negative, the overall correction to the ratio (^) in the NDR 
scheme could be as large as —20%. 

We give in Table I the values of the coefficient A(l) corresponding to a generic range 
for the form factor ratio (p0|), using two sets of Wilson coefficients C 7 Qi (m b ). The first 



set corresponds to the next-to-leading log (NLL) approximation C 7 = —0.314, Cg = 
4.154, Ciq LL = —4.261. These values correspond to the NDR scheme, and were obtained at 

H = 4.6 GeV, with rh b {m b ) = 4.4 GeV, A^ 5) = 220 MeV following Q. The perturbative 
correction from n§ multiplies CV, thus at this order consistency requires that K5 not be 
included. The one-loop matrix elements of Q\-§ are included by absorbing them into a 
scheme-independent effective Wilson coefficient Cg S defined as in [[M], 11]. At the zero recoil 



point this is given by Cf (1) = C g (m b )fj(l) +0.208 + i0.363, with 77(1) = l-0.61a,(m 6 )/7r 
0.95, where we used m c = 1.4 GeV. The corresponding numerical results for A(l) are shown 
in the first line of Table I. 

A second set of Wilson coefficients corresponds to the NNLL approximation and takes 
c nll = _ .308, C$ NLL = 4.214, C? NLL = -4.312 0, |14j. These partial results corre- 
sponding to the NDR scheme do not contain the (as yet unknown) three-loop mixing into Qg. 
The associated uncertainty in Cg NLL was estimated in [TJ]] and found to be ~ ±0.1. Also, 
the complete results for the 0(a s ) matrix elements of the Qx-q operators are not available, 
although they were recently computed in |T3|] in an expansion in powers of q 2 /m\, m^/m b . 
These approximate results are not applicable in the zero recoil region considered here. 

At NNL order the radiative correction from k 5 has to be included as well, which raises the 
issue of 75 scheme independence of the result. This has been demonstrated explicitly at NLL 
order in | JLl[j . In analogy with this result, one expects that the scheme dependence in k$ will 
be cancelled by that in the matrix elements of the four-quark operators Ya=i C^ LL (Qi) NLL 
and by that in C7. In the numerical evaluation one should use for consistency the NDR 
scheme for all quantities involved. 

The uncertainty in the value of \V ub \ extracted from ([3]) coming from A(l) is dominated 
by that in £>i(l). For |2P 1 (l)/(mj(l))| < 0.1 one finds from Table I a 3% effect in \V % 
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A precise computation of T>i(l) could help eliminate this source of uncertainty. 

In the quark model the ratio of form factors X> 1 (1)//(1) appearing in ( j30|) is always 
positive. Keeping only a S'-wave component for the vector meson wave function, one finds 
in the static limit for the h quark 

J_ (4pVb> (31) 
/(l) 6m q (40 B ) ' 1 ) 

with m q the mass of the light quark produced in the weak decay b —>■ q. The expectation 
values can be computed explicitly in the ISGW model ]^l]] with the result 



2 *i 



/(l) 2m q Pl + p 



0.093 GeV (B -> p) (32) 



x 



6 



= 0.062 GeV (B -> K*) . 

We used here the parameters of the ISGW model (3b = 0.41 GeV, m u = ma = 0.33 GeV, 
m s = 0.55 GeV, (3 P = 0.31 GeV, (3k* = 0.34 GeV. This amounts to a small positive 
contribution of 2-3% from the last term in (pOl). 





0.8 


0.9 


1.0 


1.1 


1.2 


AiVLL(l) 


-0.123 


-0.137 


-0.151 


-0.164 


-0.178 




-0.119 


-0.133 


-0.147 


-0.160 


-0.173 



Table 1. The values of the correction factor A(l) appearing in the formula (|I|) for 
a few values of the ratio (|30|) and A = 370 MeV. Two sets of Wilson coefficients 
are used, corresponding to the NLL and NNLL approximations, as explained in 
text. 

The correction to the ratio of form factors (|30"D can be also extracted from the QCD sum 



rule calculation of Ref. |2j. This S ives F( l )lfi l ) = L17 for B P and 1-18 for B -> K* . 
A similar QCD sum rule calculation in |23[ quotes the range T{\)/f{\) = 1.15^'^ for 
B — > K* . This amounts to a large positive subleading contribution of ~ 30% from the last 



term in (|3~0"D. The large discrepancy in T>\ with the dimensional analysis estimate and the 
quark model result is rather puzzling. It is not clear if this anomalously large subleading 
correction is an artifact of the sum rule computation or of the interpolation formulas in 



22j |23| . A precise determination of this formfactor is clearly important. 

In conclusion, we studied in this Letter the effect of subleading corrections to a heavy 
quark symmetry relation relevant for the determination of \V u b\ from exclusive rare and 
semileptonic B decays. The structure of the corrections to this symmetry relation is analyzed 
using the heavy quark expansion. We point out a possible large effect, depending on the 
value of an unknown matrix element of a dimension-4 operator. Lattice computations could 
eventually help to eliminate this source of uncertainty in exclusive determinations of \V u b\. 

We thank Iain Stewart for comments on the manuscript. D. P. is grateful to Andrey 
Grozin for discussions about the results of Ref. |2(J . This work has been supported by the 
DOE under Grant No. DOE-FG03-97ER40546. 
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